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A JUSTIFICATION OF THE TIMOSHENKO BEAM MODEL 
THROUGH r-CONVERGENCE 

LIOR FALACH, ROBERTO PARONI, AND PAOLO PODIO-GUIDUGLI 


Abstract. We validate the Timoshenko beam model as an approximation of 
the linear-elasticity model of a three-dimensional beam-like body. Our valida¬ 
tion is achieved within the framework of T -convergence theory, in two steps: 
firstly, we construct a suitable sequence of energy functionals; secondly, we 
show that this sequence T-converges to a functional representing the energy of 
a Timoshenko beam. 


1. Introduction 

Understanding the relation between three-dimensional elasticity and the lower¬ 
dimensional theories of elastic structures is a long-standing quest in rational con¬ 
tinuum mechanics. All classic models for shell-, plate-, and beam-like bodies rely 
on some Ansatze about their kinematic and/or static behaviour, motivated by their 
thinness or slenderness, that is, by the smallness of one or two of their dimensions. 
Such Ansatze are expedient to put together mathematical models that are both 
simple and capable to provide good-enough predictions for plenty of the intended 
applications. However, as is always the case for intuition-based models, an all- 
important experimental confirmation does not replace for a rigorous justification, 
that is, validation as a convincing approximation of an accepted parent theory. 

Justification of lower-dimensional structure theories has been attempted in a 
number of ways, some essentially analytical in nature, like the method of asymptotic 
expansion [5] and the functional analysis methods accounted for in E®, and some 
essentially mechanical, like the method of formal scaling (TO] and the method based 
on a thickness-wise expansion US]; the predictions of lower-dimensional models 
can also be corroborated by error estimates [7|- I 11 the past couple of decades, a 
noticeable amount of work has been devoted to provide a justification of various 
structure theories within the framework of T-convergence; in particular, Reissner- 
Mindlin’s theory of shearable plates has been considered in mmm- In this 
paper, for the first time, a justification via T-convergence is given for Timoshenko’s 
sheareable beam theory mi the one-dimensional theory that parallels the two- 
dimensional Reissner-Mindlin’s theory!]] 

Heuristically, justification of beam models can be achieved by calibrating the 
convergence rate of the three-dimensional elastic energy (briefly, the energy) with 
respect to the diameter of the cross section. Anzellotti et al. |T| and Bourqin 
et al. [3] were first to apply the theory of T-convergence to deduce unshear able 
beam and plate models within the setting of three-dimensional linear elasticity. 
Their results have been generalized in several ways, and different models of beams, 

J We refer the reader to [iTj for a comprehensive introduction to T-convergence, and to [2] for 
an authoritative survey of beam models. 
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plates, and shells have been studied, within the linear as well as the nonlinear 
elasticity framework. We refrain from citing the large literature on the subject; we 
only mention a paper by Mora and Muller m, where a justification of a bending- 
torsion beam model is provided within the non-linear elasticity framework. 

In the engineering community, Bernoulli-Navier’s and Timoshenko’s beam theo¬ 
ries are well accepted and constantly used in applications. As mentioned, the former 
model has been fully justified since long by means of the theory of T-convergence, 
whereas a similar justification for the latter was lacking. In our opinion, this was due 
to the fact that the “classical” procedure used in the study of dimension-reduction 
problems was excessively “constrained”. The constraints in question were relaxed 
in mm, so as to deduce the Reissner-Mindlin’s plate model; the procedure there 
used was later classified in pA| as a special case of a general scheme. Briefly, for 
a given three-dimensional problem (the so-called “real problem”), the scheme con¬ 
sists in defining a problem sequence whose variational limit approximates the “real 
problem”. 

Here is a summary of the contents to follow. In Section [2] we present the real 
problem for a linearly elastic and transversely isotropic three-dimensional beam-like 
body; in Section[3]we recall Timoshenko’s kinematical assumptions; and, in Section 
U along the lines proposed in ni and further discussed in HU, we introduce a 
sequence of variational problems, one of which is the real problem. The main 
compactness results are deduced in Section [5j while the T-limit result is stated and 
proved in Section [Gj interestingly, the T-limit turns out to be the energy functional 
of the one-dimensional Timoshenko’s beam theory. Finally, in Section [8l we show 
that the T-limit problem approximates the real problem well. 

As to notation, throughout this work Greek indices a, /3 take values in the set 
{1,2}, Latin indices i. j in the set {1,2,3}. We use L 2 (A;B) and H 1 (A; B) to 
denote, respectively, the Lebesgue and Sobolev spaces of functions defined over the 
set A and taking values in the vector space B ; in case B = R, we simply write L 2 (A) 
and H 1 (A); the corresponding norms are denoted by ||• an< ^ IHlifqA-B)- 

2. The real problem 

We consider a three-dimensional body occupying a domain under form of a right 
cylinder tt r = u r x (0, L) C R 3 of length L, whose cross-section w r C R 2 is an open 
bounded simply-connected set with Lipschitz boundary dcu r . For d r the diameter 
of the cross-section, we set e r := d r /L and we call f l r beam-like, because we take 
e r < 1. 

We denote by ei the vectors of an orthonormal basis, with vectors e a tangent to 
uj r , where the origin of the Cartesian frame we use is located, and vector e 3 in the 
direction of beam axis. For definiteness and simplicity, we stipulate that the body 
is clamped on the part <9c>fl r = cu r x {0} of its boundary. Moreover, we denote by 
U r the space of admissible displacements: 

U r = {u £ H 1 (Sl r ;R 3 ) | u = 0 on dn^r} ; 
we measure the admissible strains point-wise by means of the symmetric tensor 

Eu := i (Vn + Vu 1 ) ; 

and we denote by R 3 *^ the collection of all symmetric linear mappings of R 3 into it¬ 
self. Finally, we assume the material to be linearly elastic and transversely isotropic 
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with respect to the e 3 direction, so that the elastic-energy density per unit volume 
is given by 


^ W (E) 2 2 /i (E\i + Bf 2 ) + A (-Bn + E 22) 2 + 2T2E33 (Bn + E 22 ) 

+4/iB 2 2 + 4y (El, + Bg 2 ) + T1B33] , 

where the material moduli p., A, n, T 2 , and 7 , satisfy the following inequalities: 

(2.2) p > 0, 7 > 0, n > 0, n(A + /it) - r| > 0; 


consequently, the elastic-energy density is positively bounded below by the strain 
norm, in the sense that there is a positive constant C such that W(E) > CjB | 2 for 
every B e 

In this section and henceforth, we systematically make use of a subscript or 
superscript ‘r’ as a reminder of the fact that all kernel letters carrying that modifier 
are used in connection with a real three-dimensional equilibrium problem for a 
beam-like body. 

The elastic potential W r : lA r —>• M associated with the beam-like body f l r is: 

(2.3) W r (u) = [ W(Eu) dx, for all ueU r . 

Jflr 

Since we let fl r be subjected to a distance load b : fl r —>• K 3 and a contact load 
c r : —> M 3 , with cfyffy = dfl r \dD^r, the load potential is: 

(2.4) E r (u) := f b -udx+ f 7f -uda. 

J J d]y£l r 

Finally, the total potential is: 

(2.5) IT(tZ ) :=W r (u)-T r (u). 


Hereafter we denote by the unique minimizer of H 1 ’: 

( 2 . 6 ) i 4 in = argminn r (?l). 

U&Ar 


3. The Timoshenko Ansatz and its mechanical interpretation 

Roughly speaking, according to the kinematical Ansatz on which Timoshenko’s 
beam model is based, the beam’s cross-section is regarded as rigid, while the beam’s 
axis may deform arbitrarily. In fact, a prototypical Timoshenko displacement field 
has the following form: 

(3.1) u T = u 1 (x 3 )e 1 + u 2 (x 3 )e 2 + ( u 3 {x 3 ) + x 2 ipi (x 3 ) - x^fa)) e 3 , 

where Ui and %j) a are real-valued functions dehned on (0 ,L). 

It was shown in | 8 ] that, to within rigid global displacements, the above displace¬ 
ment field is the general solution of the following set of PDEs: 

^ ^ 2(Bw) a/ 3 — M a> /3 T U/3,a — 0, 

2 ((B'u)3q-^ ,/} — 1l 3l ap + Ua,3/3 = 0- 

The following precise kinematical interpretation of m was given: for each fixed 
axial coordinate x 3 , the first-order PDEs in (13.21) imply that cross-section fibers 
neither change their length nor change their mutual angle, while the second-order 
PDEs imply that the change in angle between an axial fiber and a cross-section 
fiber does not depend on the direction of the latter. It was also remarked in 
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that the each of the above PDEs can be interpreted as a non-dissipative internal- 
constraint condition, and the reaction stresses and hyperstresses maintaining those 
constraints were determined. 

In Q35], it was suggested that both shearable-structure theories, Reissner-Mindlin’s 
for plates and Timoshenko’s for beams, are derivable from a three-dimensional 
elastic-energy functional that includes a second-order strain contribution. This 
suggestion will be taken into account in the construction of the energy sequence 
given here below. 


4. A SEQUENCE OF VARIATIONAL PROBLEMS 

In this section we construct a sequence of variational problems, parameterized by 
e, such that the “real problem” presented in Section [2] is the element of the sequence 
whose parameter £ is equal to e r . Our construction of the typical problem in this 
sequence is achieved in a number of scaling steps, both for data and candidate 
solutions. 


4.1. Domain. For a positive, let R a := diag(a,a, 1) be a diagonal 3x3 matrix; 
moreover, for x r = (x\, £ 3 ) a typical point of f l r and for any given e > 0 , 

let and denote the sets in R 3 whose typical points are, respectively, x = 
(xi,X2,Xs) = R l t Er x r and x E = (£aq, ££ 2 , £ 3 ) = R E x (note that, consequently, 
x e = R e / £ -x r ). It follows from these definitions that 

£• 

= oj e x (0, L), where oj e := — u r , 


and that 

fl = u> x (0 ,L), where u> = u> 1 := — uj r . 

£ r 

Therefore, ui e and w are nothing but the domains in R 2 obtained by homothetic 
rescaling of co r by, respectively, factors e/e r and \/e r . Note for later use that the 
following relationships hold: 

• dv e = e 2 dv, between the volume measures of and f 2 ; 

• da E = e 2 da, between the area measures of lo e and u. 


4.2. Displacement and strain fields. Given a displacement field u e : —>• R 3 , 

we let the scaled displacement u E : id —> R 3 be defined by 

(4.1) u E (x) := R E u E (R E x); 

in Cartesian components, 

u\ = eu\ o R E , U2 — £ U2 0 R £ 5 U3 = o R E . 

It follows form definition (14.11) that 

(4.2) Vu e =R E (Vu E )R E , 

and hence that 

(4.3) Eu E = R E (Eu E )R E ] 

we call 

(4.4) 


E E u E := (i ? e ) _1 Eu E (f ? e ) _1 
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the scaled strain , and we record here its component form: 


(4.5) 


E E u e = 


/ 

(Eu e ) 11 

(Eu e ) i2 

(EU E ) is 

\ 


£ 2 

e 2 

£ 



(Eu e ) 2 i 

(■ e u e ) 22 

( E U e ) 23 



£ 2 

£ 2 

£ 


V 

(Eu e ) 3 i 
£ 

(■ e u e ) 32 
£ 

(Eu E ) 33 

j 


4.3. Elastic and load potentials, total potential. 


(i) Let 
(4.6) 


W e {E) — [2 [iE a pE a p + A ( E aa ) + 2r 2 E 33 E a 


+ 47 ^ —^ E a3 E a3 + Ti-Efg] 

(note that W Er (E) = W(E)). The elastic potential is the functional 

W E (u E ) := f W E ( E E u E ) dv 

Jn 

+ \ TR ( M W + <,3p) 2 dv, t r > 0 . 


(4.7) 


defined over the collection of all elements v £ H 1 ( Cl ; R 3 ) such that {v 3t<x p +v a ^ 3 p) £ 
L 2 (fl) (cf. definition (15.1D 1 ). 

(ii) For b : Q r —> R 3 and c r : —» R 3 , the real distance and contact loads 

introduced in Section [2j we let 


and 


"(x) = 


b r {x) := (R £r )- 1 tf(R^x) 


(R E '-)~ 1 c r (R s r x) 


on w x {L}, 

± (R Er )~ 1 c r (R Er x) on du x (0,L) 

be the scaled loads, and we assume that b r £ L 2 (fl;R 3 ) and c r £ L 2 (9/vA2; 
The load potential is the functional T : H 1 (f2; R 3 ) — \ R defined by 


(4.8) 


T 


(■ u E ) := f b r ■ u e dv + f i 

J Cl J d^Cl 


E da. 


(iii) the total potential is 

(4.9) IF ( u e ) := W e ( u E ) - T ( u e ). 

Remark 1 . We close this section by observing that energies IF and IF r coincide, 
to within a multiplicative constant. Precisely, with the use of the above data and 
solution scalings, it can be shown that 


(4.10) U r (u 6 ’') = (u Er ). 

Since the multiplicative constant e 2 does not affect the minimization process, we 
deduce that, up to a change of variables, the minimizers of IF r and IF coincide. 
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5. Compactness 


In this Section we present a compactness lemma for a sequence {w £ } with uni¬ 
formly bounded energy. The functional II E is well-defined on 
(5.1) 

A := £ H 1 (12; R 3 ) | v = 0 on w x {23 = 0} and (f3,a/3 + v a , 3/3) £ L 2 (12) } ; 

we extend II e to all of L 2 (12; R 3 ), without renaming it, as follows: 


II e (u) := 


n £ (u) for u £ A, 

00 for u £ L 2 (f 2 ; 


\A 


Moreover, for 

we let 
(5.2) 


H 1 d ( 0,L)= {veH 1 (0,L)|t>(0) = 0}, 
TV := {u £ H 1 (17; R 3 ) | 3 u?, if a £ H X D (0, L) , 


u = u° 1 {x 3 )e 1 + w°(x 3 )e 2 + (w 3 (x 3 ) +x 2 ^ 1 (^ 3 ) - xi ^ 2 (^ 3 )) e 3 } , 

We view .4 as the set of admissible displacements of a beam-like body, and TV C A 
as the subset of all displacements compatible with Timoshenko’s kinematic Ansatz. 

Lemma 1 (Compactness lemma). Let a given sequence {w £ } C A be such that 

(5.3) supll e ( u e ) < 00 . 

£ 

Then, there are a subsequence {it e } (not relabeled) and an element u of TV such 
that u e —*• u in H 1 (17; R 3 ). Moreover, 

(Eu e ) a p^-0 otL 2 (17). 

Proof. By means of Holder’s inequality, (12. 2D . the trace theorem for //^-functions, 
and Korn’s inequality, which holds because u E = 0 on ui x {x 3 = 0}, we find: 

n Hu e )>CJ 2 [ ((^u E )lp + e^E E unl 3 + (E E u E )l 3 + ^{ul a p + u E a>3 p) 2 )dx 

a,f> JnX - £ J 


-m 

> Cl II Eu 11^2(0) — C2 IIif 1 (O) 
— ^3 ||W lljyi(n) ~ C 2 ll' u ' e II(O) > 


l 2 (n) ll u llL 2 (n) \\ c 


li 2 (a D n) 


llL 2 (a D n) 


> C4 ||w e || ff i( 0 ) C5. 


It follows from assumption (15.311 that {w £ } is a bounded sequence in H 1 (12; R 3 ). 
Moreover, from the first line of the inequalities above we deduce that all sequences 


{Eu e ) 


a/3 


> {( Eu l a3 }> {(^ £ ) 33 }, 


U la0 + U a ,3p 


are bounded in L 2 (17). Hence, up to a subsequence, u e ^ u in H 1 (17;R 3 ) for 
some u £ H 1 (17; R 3 ). In addition, from the fact that | < ‘ Eu J ap j is a bounded 
sequence in L 2 (17), it follows that u e a ^ + Up a = 2 ( Eu e ) a p -A 0 in L 2 (17); and, 
u a 0 + a u ct,P + u p,a in L 2 (17); we conclude that u a< p + up, a = 0. Finally, 
with a similar argument, we can show that u 3 a/3 + u e a 3 p — 1 u 3 , a p + u a , 3 p = 0 in 
L 2 (17). Thence, u £ TV. □ 
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6. The T-limit of {IT} 

In this section we analyze the T-limit of the sequence of functionals IT. For 
Wf {E i3 ) := min { W e ( V' E a3 (e a ® e 3 + e 3 < 8 > e„) + E 33 e 3 <g> e 3 

9c.fi 

ot 

i \ ’ 9a/3 r i \ l 

+ 2^ 9 l e « ® e /3 + ep < 8 > e a \)), 

a,yS 

a simple calculation shows that 

(6.2) Wf (E l3 ) = \[±l{-)\El 3 + El 3 ) + (n - yt-) £ 3 3 ]' 

2 £ r A + fi 

Define 

(6.3) W T (S i3 ) := W}(E a ) = 1 -[^e~ 2 {E 1 3 + E 2 23 ) + ( n - ^)£&], 
so that 


(6.4) IT T e (T i3 ) = Wr(eE 13 , eE 23 , T 33 ). 

Define W r : L 2 (12; R 3 ) —> R by 


(6.5) 




fa w r (( Eu ) i3 ) dx > 

00 , 


u G TT 
otherwise 


and II : i 2 (12;R 3 ) -> R by 


( 6 . 6 ) 


n(u) 


w r (ti)-j(«), «en> 

oo, otherwise 


Theorem 2. The sequence {II e } T-converges to II in the L 2 (12; R 3 ) -topology, that 
is to say, 

( i ) (liminf inequality) for every sequence {u e } C L 2 (12; R 3 ) and for every u G 
L 2 (12; R 3 ) such that u e u in L 2 (12; R 3 ), 

(6.7) hminf IT (u £ ) > n(u); 

£->0 

(ii) (recovery sequence) for every u G L 2 (12; R 3 ), there is a sequence {u e } C 
L 2 (12; R 3 ) such that u s -A u in L 2 (12; R 3 ) and 

(6.8) limsuplT(T) < n(u). 

£—>-0 

Proof. We start by proving (i). Let {u E } C L 2 (12; R 3 ) and u G L 2 (12; R 3 ) such 
that u e -A- u in L 2 (12; R 3 ). The inequality in (16.71) is nontrivial only for 

liminf IT (u e ) < oo. 

E-m 

Thus, (by passing, if needed, to a subsequence) we may assume that 

lim IT ( u e ) = liminf IT ( u E ) < oo. 

£—>0 £—>-0 
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By Lemma [U the sequence {it £ } converges weakly in H 1 (52; R 3 ) to an element 
u £ TV. Thus, 

liminf II £ (u £ ) > liminf { [ W E (E E u E ) dx - T (u E )}, 

e^O e^O J n 

> liminf {[ W E ((E E u E ) i3 )dx-T(u E )}, 

= liminf{ f W T (£ (E E u E ) 13 ,£ (E E u E ) 23 , (E e u E ) 33 ) dx - T (u E )}, 
£->-u Jn 

= lim inf { f W r ((Eu E ) 13 , (Eu E ) 23 , (Eu E ) 33 ) dx - T (u e )}, 

> [ W T ((Eu) 13 ,(Eu) 23 ,(Eu) 33 )dx-T(u), 

J n 

= W T (u) - T(u) =n(u). 


We point out that: in the first inequality we dispensed with second-order terms; 
the second inequality makes use of (16.11) ; in the third line we applied (16.41) ; the last 
inequality follows from the convexity of W T and the continuity of T. 

We now prove (ii). Let u £ L 2 (52; R 3 ). Since inequality (16.81) is trivial for u 
TV , we only need consider the case when u £ TV, and hence has the representation 
specified in (15.21) . 

At first, we restrict attention to functions u 3 , ipa belonging to C°°(0,L) and 
equal to zero in a neighborhood of x 3 = 0 ; note that these functions form a dense 
subset of Hp (0, L ). We consider the sequence whose typical term is: 

u E := u + e 2 u, 


where u is defined by 

ui := 

(6.9) u 2 := 

U3 ■■= 


X ^ 

+ x l x 2^1,3 - ^-^2,3 + yV^), 
X 2 X 2 

-t/(x 2 u 3 3 - X 2 X 1 1p 2 ,3 + y V>1,3 - yt/’i.s), 
0 , 


with r) := t 2 I 2 (/j, + A). A simple computation yields: 


Thus, 

W £ ( u E ) 


(-E E u E ) n 

= (Eu) n = -77 (Eu) 33 

(E e u e ) 12 

= 0 , 

(E e u e ) 22 

= (Eu) 22 = -V (Eu) 33 

(E E u E ) a 3 

= (Eu) a3 + £ 2 (Eu) a3 , 

( E u ) 33 

= (Eu) 33 - 

j 2 ( Eu ) 33 

+ 4A t) 2 (Eu) 33 T 2 t 2 (Eu) 


, , 2 , £ . 2 1 / (Ell) a3 + £ 

+n (Eu) 2 3 + (—) 2 47 (--^4- 

Cr 


2 {Eu) a 3 2 


+ rX Er c £ ) 2 g 4 (^3 ,a.0 + Ma,3/3) 2 ] dx. 
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Upon rearranging, we obtain: 

7-2 


w £ 


<“•>=/> 


4t) W 33 + (^) 2 47 ((£«) «3 + £ 2 (£fi) q3 ) 2 
/l “T A £ r 

H - ^ ^ (^r £ (^3,a:/3 H“ ^o;,3/?) ] dx. 

a,(3 


By passing to the limit, it follows that 


limVP(P) = [ h(l )H 1 (Eu)l 3 + (T 1 -^-)(Eu)l 3 ]dx = W r (uy, 
e->0 in 1 £ r /i + A 

since 

lim T (u e ) = T(u), 

E-S-0 

the proof of ( ii ) is achieved, under the smoothness assumptions stated at the be¬ 
ginning of this paragraph. 

The general u G TV case is handled by a standard diagonalization argument. 
Indeed, let {life} C TV n C°° (U;R 3 ) such that Uk — > it in H 1 (Q;R 3 ); moreover, 
let {u|} be the recovery sequence for Uk as defined by (16.91) . Since 

lim lim |K -u||tfi ( n) = 0 , 

fc-J-OO E—>0 


and 


lim lim IP (u E k ) = lim II (uk) = II (it) 

k—>oo e —>0 k—too 


we can find an increasing map e —> k e such that —>• u in H l (Ll] R 3 ) and 

lim sup IP (ui ) < lim lim IP (u E k ) = II (it) . 

£ ^.q e k — y OO £—^0 


□ 


Theorem 3. Let u^ri n be the minimizer o/IP, and let it m i n be the minimizer of II. 
Then, 

(i ) the sequence IP (itmin) converges to II (u m in); 

(ii) the sequence {itmin} converges to it m i n strongly in H 1 (fl;R 3 ). 

Proof. As to the first claim, given that sup £ {IP(«P n )} < oo, it follows from Lemma 
[□that {Uminl converge weakly in H 1 (fl;R 3 ) (up to a subsequence) to iti; m G TV. 
For a general u G L 2 (U; R 3 ), by (ii) of Theorem^ there exists a recovery sequence 
{P} such that 

H(u) > lim SUp IP (u E ) > lim SUp IP (itmin) • 

£—>-0 £—^0 

From (i) of Theorem [2j we deduce that 

(6.10) n(zt) > lim sup IP (itmin) > lim inf IP (uP n ) — ^ (itiim) ■ 


£—>-0 


£—>-0 


Now, u can be chosen arbitrarily; we take u = u m i n , and obtain II ( u m [ n ) > II (wii m ); 
as n (itmin) < n (it lim ), part ( i) of the theorem follows and itmi n — L it m i n in 
H l (U;R 3 ). 
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As to the second claim, by means of m , <ED, and (jG.dj) , wg deduce that 

W" (Umin) — ( M min) > / [IK E (E E M^in) — W"r ((-E/Wmin)^)] dx, 

J O 

( 6 . 11 ) = [ [W e (E s u £ min )-W E {(E E u min ) i3 )\dx, 

Jn 

> [ [W £ ((E E u E min ) i3 )-W E ((E E u min ) i3 )}dx. 

Jn 

On recalling the form of W E given in (16.21) . it follows after some computations, with 
the aid of ( 12 . 21 ) . that there exists a constant C > 0 for which 

[ ((E E u E min ) i3 ) - W E ((E E u min ) i3 )} dx, 

Jn 

— ( y ' ^11 (K'' u, min)i 3 — (Etlmin) j 3 || £2 + J [(Eu m - ln ) i3 — (Eu m in) i3 ] {ElLmin) i3 dx'j . 
Combination of this inequality with (16.111) yields: 

11" ( u min) — n («min) ^ C ^|| (El«mi n ) i3 ~ {Eu m i n ) i3 \\ L 2 

i 

+ J [(EM m i n ) i3 — (Ezi m in)^ 3 ] (Eu m i n ) i3 dx'j 

+ JF (Umin) F (llrnin) ■ 

As u e m - m w m i n in H 1 (fl;R 3 ), we deduce that 

/ [(Eu m i n )i 3 — (EM m in) i3 ] (£u m in) i3 dx 0, 

Jn 

thence, by part (i) of the theorem, we find 

0 > limsup^] \\(Eu £ min ) i3 - (Eu min ) i3 \\\ 2 . 


£—^0 


Thus, {Eu^^) i3 —>• (Eu m i n ) i3 in L 2 (fi). As the sequence {II e (w^in)} i s bounded, 
it follows from Lemma □ that {Eu e min ) a p -A 0 = {Eujn in ) a/3 in L 2 (fl); hence, 
Eu^ in —>• Eumin in L 2 (fl;K 3x3 ). An application of Korn’s inequality concludes 
the proof. □ 


7. The T-limit potential in terms of limit displacements 

Once again, recall that, for every u £ TV, there are u°, i/) a £ Hp (0 ,L) such 
that 

(7.1) u = Ui(x 3 )ei + U 2 (x 3 )e 2 + (w 3 (:r 3 ) + x 2 tpi(x 3 ) - xiip 2 (x 3 )) e 3 
(cf. (15.21) 1. It follows that the non-null associated strain components are: 

(Eu) 13 = ^{u° 3 -1p 2 ), (Eu) 23 = i(U2,3 + V’l) ! {Eu) 3 3 = U 3t3 +X2lpl,3-X 1 ll>2,3- 
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Thus, in view of (16.311 . the elastic potential (16.51) reads: 

dx 3 j (7£r 2 [(Wl.3-V’2) 2 + (M2 i 3+V’l) 2 


+(ri - 


■)( m 3,3 + a ; 2'01,3 - Xi1p2,3) 2 ^jda , 


A + n 

whence, on choosing for the first two Cartesian axes the principal axes of inertia of 
the cross-section w, we arrive at the one-dimensional elastic-energy functional of 
Timoshenko’s beam theory: 

r L 

W T (u) = ' 


(7.2) 


where 


-0) = \ [ (^y[K,3 - V>2) 2 + («2,3 + ^i) 2 ] 

z JO K £ r 


+ (n - 


A + /j, 


)[A(lt 3 j3 ) + + </ 2 ' 02 j 3 ]^ dx 3 


A := 


:= / da, J\ := / x 2 da, J 2 := / 

J UJ j UJ J UJ 


x? da. 


denote, respectively, the area and the moments of inertia with respect to the Xi 
and X 2 axesQ 

It also follows from (E3D that the load potential (14.81) takes the form: 


F (u) = f ( b\u ° + b 2 u 2 + b 3 u 3 + x 2 b 3 ijJi — Xib 3 ip 2 ) dv 
Jn 


L 


(c>? + c 2 u 2 + C 3 M 3 + X 2 C 3 V 4 - XiC 3 ^> 2 ) da, 

I dj^Q, 

whence we deduce the one-dimensional load functional of Timoshenko’s beam the¬ 
ory: 

r L 

F(u) = 


(7.3) 


( u ) = / (hu° + f 2 u° + /3W3 + "iiV’i + ^2) dx 3 
J 0 

+ Fiu°(L) + F 2 u° 2 (L) + F 3 u° 3 (L) + AhML) + M 2 MQ, 
where, for almost every X3 £ ( 0 , L ), we have set 

fi(x 3):= / &C(xi,x 2 ,x 3 )da+ / c[(xi,x 2 ,x 3 )ds, 

•J uj J duj 

mi(x 3 ) := / x 2 b 3 (xi,X2, x 3 ) da + / x 2 c 3 (xi, x 2 ,x 3 ) ds, 

•J uj J duj 

m 2 {x 3 ) := — / xi6 3 (xi,x 2 ,x 3 ) da — / XiC 3 (xi, x 2 , x 3 ) ds. 

J uj J duj 

and 

.Fi := / c r i {x 1 ,X2,L) da, Mi := / x 2 c 3 (xi, x 2 , L) da, M 2 := — / xic 3 (xi, x 2 , T) da. 

J UJ j UJ j UJ 


2 Under the assumptions, 


/ 

./ a; 


( da = 0 and 


/ 

•/ Co> 


011 X 2 da = 0. 
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Interestingly, for the elastic potential as in (17.21) and the load potential as in 
(17.31) , the total potential (16.61) may be decomposed into an axial-stretching part II a 
and a bending part Ilf,: 

(7.4) n(u) = n a (u§) + n 6 (u“,^a), 
where 

(7.5) n a (ug) := i(n - —j— )A [ « 3 ) 2 dx 3 - [ f 3 u° 3 dx 3 + F 3 u° 3 {L ), 

z A + p j 0 j 0 

and 

n 6 (w“,W:=^ [ (u ° lt3 - ^ 2 ) 2 + + ip 3 ) 2 dx 3 

z e r J 0 

( 7 - 6 ) + \{tx - y~~) f - 71 ^ 1,3 + J ^l, 3 dx 3 

Z A + Jt Jo 

flU°i + / 2 U 2 + TO 1^1 + TO 21 A 2 rfa; 3 
- Fiu°(L) - F 2 m°(L) - Miipi(L) - M 2 ML)- 

Thus, the minimization problem splits into two independent problems: the one for 
m 3 , to be determined by minimizing the axial-stretching potential n o , the other for 
ur a and ij) a , to be determined by minimizing the bending potential Ilf,. 

Remark 2. In the isotropic case, the elastic-energy density defined in (12.11) only 
depends on the two parameters A and /z, because we have that 

7 = M, Ti = A + 2/z, r 2 = A. 

Consequently, the elastic modulus appearing in both the axial-streching and the 
bending potentials reduces to the Young modulus of the material: 

r| 3 A + 2 fi) 

T l — - = -. 

A + A + )jl 

Remark 3. The parameter e r enters both the three-dimensional limit elastic- 
energy density (16.31) and the elastic part of the one-dimensional bending potential 
(17.61) (but not the elastic part of the axial-stretching potential (17.51) 1. If, ceteris 
paribus, we were to let e r —> 0 in (16.31) . we would achieve effortlessly a justification 
of the Bernoulli-Navier beam model, whose total potential obtains by letting e r —> 0 
in (17.61) . In fact, in the envisaged limit, the shear strains E a 3 would be forced to 
converge to zero, which is tantamount to take 1 /I 2 = «i 3 and %/ji = — 3 in dm 

This remark supports the engineer idea that the Bernoulli-Navier model is fine for 
very slender beams, whereas the Timoshenko model is preferable whenever beams 
are not so slender. 



8. Summary and conclusions 

In Section [2j we have defined the total potential IT of a linearly elastic three- 
dimensional beam-like body; we have denoted the minimizer of this potential by 

<4in = arg min IT(m); 

UEU r 
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and we have denoted the ratio between the diameter of the cross-section and the 
length of the beam by e r , a parameter that measures the slenderness of the beam¬ 
like body at hand. In Section 3, we have recorded the form of the Timoshenko 
displacement field, and we have briefly discussed its mechanical significance. In 
Section [4] we have constructed an e-sequence of functionals IP, such that IP r is 
proportional to the “real” functional IP (see (14.101) 1. After studying, in Section 0 
the compactness properties of sequences with bounded energy, we have identified 
the T-limit of the sequence IP in Section [HI The T-limit II turns out to be the 
total potential of a Timoshenko beam; in Section 7, we have written it in terms of 
the fields that parameterize the class of Timoshenko’s displacements (see (16.Gl) and 
(17.41) - (17.61) 1. and we have shown that these parameter fields can be determined by 
solving two independent minimum problems for, respectively, axial stretching and 
bending. 

In accordance with the notation introduced in Theorem [3J let 

u min = arg min IP (u) and u m i n = arg min II(m). 
uEA u£TT> 

The minimizer u m - m is a Timoshenko-type displacement. Since II r essentially co- 
incides with IP r , we deduce that u( nirl essentially coincides with u^ in . The exact 
relation between them follows immediately from (14.101) and (14.11) . and is 

(8.1) Uminfc) = Rer u r min (R er x). 

In Theorem [3j we have shown that 

u E min ->• itmin in U 1 ^;® 3 ); 

thus, we can loosely say that, for small e, u E n - m is well approximated by it m i n . In 
particular, for e r very small, it^ in is well approximated by u m i n ; we concisely write 
this as u(P n « u min . By m we therefore find an approximation of the “real” 
displacement u( nin , namely, 

uJninOO ~ (R er )~ 1 u ln i n (R 1 /Sr x r ). 

This relation states that the “real” displacement in is well approximated by a 
Timoshenko-type displacement and that such an approximation can be constructed 
with the use of the minimizer of the T-limit we found. 
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